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Abstract 

Once chosen the dynamics in one frame, the rest frame in this paper, the Bakamjian and Thomas method 
allows to define relativistic quark models in any frame. These models have been shown to provide, in the infinite 
quark mass limit, fully covariant current form factors as matrix elements of the quark current operator. In this 
paper we use the rest frame dynamics fitted from the meson spectrum by various authors, already shown to provide 
a reasonable value for p 1 2 . From the general formulae for the scaling invariant form factors ^ n \w), and 

Tyl(w), we predict quantitavely the B semileptonic branching ratios to the ground state and orbitally excited 
charmed mesons D, D* and D**. We check Bjorken’s sum rule and discuss the respective contributions to it. We 
find f(w) ~ (2/(1 + w)) 2 , resulting from the fact that the ground state wave function is Coulomb-like. We also find 
T3/2 — 0.5(2/(l+u>)) 3 and ti/ 2 (w) -C T3/2 (w)- Very small branching ratios into j = 1/2 orbitally excited D's results. 

The overall agreement with experiment is rather good within the present accuracy which is poor for the orbitally 
excited charmed mesons. We predict a ratio Br(B —> D^lv)/Br(B —> D\lu) = 1.55 ± 0.15 as a mere consequence 
of the heavy quark symmetry. If some faint experimental indications that Br(B —> Dilu) ~ Br(B —> D^lv) were 
confirmed, it would indicate a sizeable 0(l/m c ) correction. 
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1 Introduction 


The experimental progress in beauty physics these last years is really astounding. Subdominant decay channels such 
as the decay into non charmed final state are now currently seen and measured with an increasing accuracy. In this 
paper, we are focusing our attention on another class of non dominant decays, namely B —> D**. Their suppression 
relative to the dominant B —> is of course not due to CKM angles but to the dynamics of the decay, and mainly 

to the orthogonality of the orbitally excited spatial wave function of the D**’’ s to B' s the ground state one. 

Among the four non-strange expected D**' s, two have been seen, the £)|(2460) ± with I(J P ) = |(2 + ), and the 
Hi(2420) ± with I(J P ) = ^(1 + ), which have rather small widths (~ 20 MeV). This small width is usually understood 
as a consequence of being j = 3/2 states, j being the total angular momentum of the light quanta (light quarks and 
gluons) of the system! 2 3 !. The small widths of these states is attributed to the centrifugal barrier effect on their decay. 
Indeed, the dominant decay channels of the D**’s are and D**(j = 3/2) —> D^*\j = 1/2)7r is only allowed 

through a L = 2 partial wave by parity and j conservation. Of course j conservation is only valid at leading order 
in 1 /m c . It appears that the main expectations of HQET seem, up to now, to be rather well fulfilled, even though 
the charm is not such a heavy quark and important corrections would not be a surprise. According to what seems 
to be the standard use, we will denote the j = 1/2 states: Z)*(~ 2420) with I(J P ) = ^(1 + ) and Dq(~ 2360) with 
I{J P ) = \{ 0+). 

There are altogether three main reasons which lead us to study these D**' s with some care. First, narrow 
resonances are in general preferred for phenomenology for many practical reasons: they are better isolated from the 
background and from other resonances, the very convenient and frequent “narrow resonance approximation” can be 
rightfully used, etc. Second, the apparent success of HQET at dominant order in l/m c leads us to hope a reasonable 
success when describing the B —> D**lv decay in the infinite mass limit without having to face the large number of new 
operators and unknown parameters at next order in l/m c . Last, but not least, the excited hadrons represent a domain 
in which, to our knowledge, quark models remain the only tool to make predictions, beyond general consequences from 
Heavy Quark Symmetry (HQS). Lattice QCD as well as QCD sum rules meet huge difficulties to deal with hadrons 
beyond the ground states. 

Thus, B —> D** is the privileged domain to test a class of relativistic quark models a la Bakamjian and Thomas 
(BT) ||], which we have recently proven Q] to provide a covariant description of semileptonic decays in the heavy mass 
limit. It was shown (d) that this class of models satisfies fully HQS and that the p 2 Isgur-Wise slope parameter was 
bounded in this class of models: p 2 > 0.75. It also satisfies, [||, 0 , duality properties amounting in particular to the 
Bjorken sum rule 0]. In | jTo| it is shown that they satisfy the new sum rules involving annihilation constants that 
we have proposed in |j)|. 

The method of BT proceeds as follows. Given the wave function in, say, the rest framej], the hadron wave 
functions are defined in any frame through a unitary transformation, in such a way that Poincare algebra is satisfied. 
The mass operator, i.e. the Hamiltonian operator in the rest frame, is not specified in the BT approach. The mass 
operator is simply assumed to be rotationally invariant, to depend only upon the internal variables (to commute with 
total momentum and center of mass position) and, of course, to conserve parity [Q. In other words, the spectrum of 
heavy flavored mesons is left free, and all the above mentioned properties are valid for any spectrum. 

These welcome properties make the BT method appealing enough to try real phenomenology. The hadron 
spectrum has been studied many times, and several good fits to the masses of heavy flavor mesons exist in the literature. 
The question we ask is, given any mass operator which fits well the mass spectrum, what does it predict for the slope 
—p 2 of the Isgur Wise function at v.v’ = 1, and for the B —>■ D**lv decay ? The results about p 2 have been published 
and discussed in a separate publication [|j] . It has been stressed that the detailed shape of the potential (more or less 
singular at the origin, etc.) did not influence much p 2 1 while the form of the kinetic energy (KE) plus mass, whether it 
is relativistic (\Jp 2 + m 2 ) or non relativistic (to + p 2 /2m), did change very sensibly the result. While the relativistic 
KE gives p 2 = 1.00 ± 0.05 for all models that we have considered, the non relativistic KE provides much larger values, 
p 2 ~ 1.3. Theoretical and experimental estimates of p 2 clearly favor the value close to 1, i.e. the relativistic KE 0 . It 
was a good surprise for us that the more sensible relativistic KE combined with a covariant way of boosting the states, 
the BT method, provides a relatively model independent and quite reasonable value for p 2 . It was a good surprise 
because the relatively stringent lower bound p 2 > 3/4 [Q, valid for any mass operator, could induce the fear that a 
spectroscopically successful mass operator, and its eigenfunctions which necessarily are very different from the ones 
saturating the lower bound, would give a much too large p 2 . This fear could be reinforced by the fact that gaussian 

2 If the DJ has to be a pure j = 3/2, the two expected J = 1 states with j = 1/2, 3/2 should mix due to 1 /m c corrections to the Heavy 
Quark Effective Theory (HQET). This mixing is estimated to be small jlj and the observed state to be dominantly j = 3/2. 

3 Any starting frame can be chosen, for example the infinite momentum frame although some care is then needed. 
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wave functions give a large value p 2 > 1.2. It was quite gratifying to realize that the properties of the spectrum and 
the short distance potential predicted from perturbative QCD, plus the relativistic KE was shaping the wave function 
towards a lower and better value of p 2 , and towards a larger value of the wave function at the origin. The latter effect 
and the related leptonic decay constants are discussed in jnjj . 

To summarize, we have a model which has several welcome properties, and provides a very reasonable value 
of p 2 , and which, among its beauties, satisfies the Bjorken sum rule Q, Q which, in particular, relates the value of p 2 
to the functions which enter in the B —> D**lv decay, eq. (4.1). 

This is not a trivial achievement. It comes in this class of models because the boost of the wave functions 
is a unitary transformation that keeps the closure property of the Hilbert space in all frames. It is then more than 
tempting to make predictions with this model for B —> D**lv precisely at a time when these branching ratios start 
being measured. 

In j(J we have derived the general formulae giving the HQET form factors t \ i 2 { w ) and t 3 / 2 ( w ), which are 
relevant for B —> D**lv, from the eigenfunctions of the mass operator. 

The content of the BT method has been explained in some detail in |Z| and (Q. It will not be repeated here. 
In this paper we will use the same set of mass operators fitted to the experimental meson spectrum which has proven 
to provide a reasonable value of p 2 |ll|]. In section 2 we describe the four models from which we borrow the mass 
operator, and our diagonalisation procedure. In section 3 we derive the formulae for the partial widths. In section 4 
we check the Bjorken sum rule. In section 5 we describe the shape of the functions f;(w) and Tj(w). In section 6 we 
will present our predictions for semileptonic branching ratios and compare them with experiment and other models. 
Finally we conclude. 


2 Diagonalization of the mass operator 

2.1 Mass operators 

The method of Bakamjian-Thomas || [|j provides a fully covariant J|] description of the current matrix elements in 
the infinite mass limit, once given a mass operator, i.e. a Hamiltonian operator in the meson rest frame. For the 
latter the only condition that we impose is that it has to describe correctly the meson spectrum. We did not try to 
invent our own Hamiltonian but rather resorted to literature. We have chosen four quark models: Godfrey and Isgur 
(GI) [ p~2| , Veseli and Dunietz (VD) |l3|], Cea, Colangelo, Cosmai and Nardulli (CCCN) Q and finally Isgur, Scora, 
Grinstein and Wise (ISGW) |lGj]. The GI model is rather complex and because of that, one may not agree on all the 
ingredients or assumptions which enter it, or one may find it difficult to evaluate their respective impact on the final 
result. Nevertheless, what makes the model outstanding is its covering of the whole spectroscopy, from light to heavy 
quarks. The GI model includes relativistic features, among which the square root kinetic energy: 


K = \Jp 2 + m\ + \Jp 2 2 + ml , ( 2 . 1 ) 

Starting from the linear plus Coulomb potential complemented with spin-spin and spin-orbit forces, GI apply a 
complicated procedure to smooth the singular parts of the potential (using for example an a s (Q 2 ) which is finite when 
Q 2 —> 0) and smearing the potential via momentum dependent terms meant to mimic relativistic effects. Since we 
consider in this paper the infinite mass limit m\ —> oo, we have thrown away in the present calculations all terms 
which vanish in that limit, except when checking our numerical code against their spectrum. We use the following set 
of parameters in the notations of |b| : 

b = 0.18GeV 2 , c = -253MeV, cr 0 = 1.8GeV, s = 1.55, e SO (V ) = -.035, e SO (S) = .055 (2.2) 


a s (Q 2 ) = 0.25e-« 2 + 0.15e-« 2 * / 10 + 0.20e-« 2 / 1000 , 7 r = 0.5, 72 = 73 = 


m-i = 10 4 GeV, m 2 = 0.22 GeV. 


where the mass mi = 10 4 GeV is taken as a good approximation of TOi = oo Q 

The VD model is a genuine linear plus coulomb model fitted to the beauty and charmed flavored mesons: 


4 Q!q 

H = K -b br + c 

3 r 


(2.3) 


4 When computing the B meson decay widths, we will use the amplitudes computed at m\ = 10 4 GeV but the phase space factors will 

be taken with the physical masses of the B , the and D **. 
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where K is the relativistic kinetic term (2.1) and we use 


a s = 0.498, b = 0.142 GeV 2 , c = -350MeV, mi = 10 4 GeV, m 2 = 300MeV 
The CCCN model (13] uses a potential a la Richardson p5[| 


H = K 


8tt 


with K in © and 


4 

m = - / 

77 Jo 


33-2 n f 

sin(gf) 


A Ar- 


/(Ar) 
A r 


dq 


1 


1 


q VMl + 'l 2 ) 


(2.4) 


(2.5) 


( 2 . 6 ) 


Richardson’s potential has the property to be linear for large r and to behave as predicted by asymptotic freedom for 
r —y 0: 


V(r)- 


8tT 


1 


33 — 2 nf r ln(Ar) 


(2.7) 


We use the set of parameters: 


A = 0.397 MeV, n f = 3, mi = 10 4 GeV, m 2 = 38 MeV (2.8) 

and, as in [fj~i| , we truncate the potential to a constant value, 

V(r) = V{r m ), forr < r m , with r m = ^^GeV -1 (2.9) 

5.33 

Finally, the ISGW model Jh|, contrary to the three previous ones, has a non relativistic kinetic energy term 


H = 



_P2_ 

2 m 2 


4a s 
3 r 


br 


( 2 . 10 ) 


with the set of parameters 

a s = 0.5, b = 0.18GeV 2 , c = -840MeV, mi = 10 4 GeV, m 2 = 330MeV (2.11) 


In reference ||l6| form factors are computed according to the non-relativistic dynamics, with some adjustment 
of the parameters. We insist that in this paper we do not use the latter estimate, but we only take from [^6| the 
Hamiltonian at rest and compute the form factors according to BT method. The result is indeed totally different. 
The label ISGW refers to BT estimate of form factors from ISGW Hamiltonian, while, when necessary, we will refer 
to ISGW’s computation of form factors as ISGW-nonrel. 

2.2 Diagonalization 

For each of the four mass operators we have diagonalized the sectors l = 0 and l = 1 using a basis of harmonic 
oscillator eigenfunctions which is truncated by keeping only the n max + 1 lowest states. The code first computes 
the matrix elements of the Hamiltonian in this truncated Hilbert space, and then diagonalizes the matrix keeping 
the eigenvectors and the eigenvalues. We have performed all the calculations with n max = 10,15, 20 and for the GI 
model with n max = 24. Before sending the mass of the heavy quark to 10 4 GeV, we have checked, using the masses 
mb,m c fitted by the authors, that we had good agreement with them for the meson masses. In the case of GI the 
check needed the introduction of spin dependent terms which disappear in the infinite mass limit. VD has singular 
eigenfunctions for r —> 0 due to the coulombic singularity. CCCN cut the coulomb-logarithmic singularity at r m , but 
a discontinuity results in the potential at r m . Therefore, VD’s eigenvalues in the l = 0 sector converge rather slowly 
when n m ax —> oo, and this is also the case to a smaller extent for CCCN. This does not prevent a rapid convergence 
of the matrix elements which we will consider in this paper. Only the leptonic decay constants are infinite and the 
sum rules of || are ill-defined for the VD potential [[To| . 
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3 Analytic formulae for B —> D ** decay 

3.1 €(w), T 1/2 (w) and t 3/2 (w) 

It is well known, that in the heavy quark limit the current matrix elements between B and the even parity D** 
states (odd parity D W) can be expressed as functions of two (one) universal functions, Tj(w) (£(w)), where w = v ■ v' 
and where j = 1/2, 3/2 is the total angular momentum carried by the light quantas (light quarks and gluons) in the 
final states, j is a conserved quantum number in the infinite mass limit. In Q (jp), we have derived the expressions 
for t 1/2 (w) and t 3/2 (w) (£(«;)) in the BT type of models, once given the eigenfunctions of the mass operator: 


^ ^ 1 /' dp 2 \J(p2-v')(p 2 -v) P2-W + v) + m 2 (v'.v + 1 ) , / 2 2 * 

= l(i?—3- y( K v + mg )( K ., + TO ) y((K - ,,) v(fe ' 


v) 2 — m 2 ) (3.1) 


t 1/2 (w) = 


1 f dp 2 1 \/{p 2 ■ v')(p 2 ■ v) 


2^3 


and 


t 3 / 2 {w) = 


<t>l((p 2 -v) -m 2 ) p{{p 2 .v) to 2 ) 


(2tt) 3 p° ^(p 2 . v > + m 2 )(p2 ■ v + m 2 ) 

(p 2 ■ v)(p 2 ■ v' + m 2 ) — {p 2 ■ v')(p 2 ■ v' + (v ■ v')m 2 ) + (1 — v ■ v')m\ 


dp 2 1 \/{p2 ■ V')(jP2 ■ v ) 


\/3 1 - (v ■ v') 2 J (27r) 3 p 2 (j> 2 ■ v' + m 2 ){p 2 ■ v + m 2 ) 


1 — v • v' 


^^ ( (P 2 ) -m 2 ) ^((P2-w) -m 2 ) 


1 


2 1 + v • v 


7 [P 2 • (v + v')\ - {p 2 • v)(2p 2 ■ v' - m 2 ) - (p 2 ■ v') [p 2 • v' + (v ■ v')m 2 ] - 


1 — V ■ v' 


(3.2) 


(3.3) 


where tp(p 2 ) (4>j(p 2 )) is the radial part of the wave function of the B (D**’s). A detailed account of our conventions 
can be found in (fll. 


3.2 Differential decay widths 

Leaving aside radial excitations, there are four l = 1 states, the D**, that we will label according to the common use: 
Dq, D\ with j = 1/2 and D i, D 2 with j = 3/2. The lower index stands for the total angular momentum J. Two 
experimental resonances, the Z?i(2420) and the D 2 (2460), are considered as good candidates for the j = 3/2 states. 
In the infinite mass limit, to which we stick here, the mixing between both J = 1 states (j = 1/2 and j = 3/2) is 
forbidden. The 0(l/m c ) mixing seems to be small [jp. 

To make short a long story, we have checked the formulae in jp] with which we agree except for three misprints: 
in eq.(45) in jp one should read in the denominator 327r 3 instead of 327 t 2 and r = rnx Qj /mp Q . (for example r = 
Too** /ms) instead of r = m ‘x Q / m ^ Q ■ Finally we agree with table III in jp if one defines 5 = [2 rw 2 — (1 + 4r + 
r 2 )w + 2(1 + r 2 )]/4r instead of <5 = [2 rw 2 — (1 + 4r + r 2 )w + 2(1 + r 2 )]/2r. 

While we have computed the wave functions and the universal functions Tj(w) in the infinite mass limit, the 
kinematics is taken with the physical masses of the particles (or, in the absence of measurement, some estimated 
physical masses), i.e. 


ms = 5.279, mp = 1.865, tod* = 2.010, mo* = 2.360, mu* = 2.420, mp 1 = 2.422, = 2.459, (3.4) 

We also use 

V cb = 0.041, G f = 1.16638910~ 5 GeV. (3.5) 

and r = mx/friB , X being the charmed final meson. The calculation of the differential decay widths is then standard, 
leading to (taking the mass of the final lepton to zero) 


^(B^Dlv) = I ^6| 2 ^^r>(t«, + l)(^-l) 1 / a (ti,-l)(l + r) 9 |f(t«,)| a 


(3.6) 
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dr 

dw 

l B- 

-> D*lv) = 

= \v cb \ 2 

G 2 F m% 

487T 3 

r 3 (w + 1 )(w 2 

dr 

dw 

l B- 

-Dllv) -- 

= \Vcb\ 2 

G 2 F m 5 B 

487T 3 

4 r 3 (w 

- !)(^ 

dr 

dw 

C B- 

-> D\lv) = 

= IKfcl 2 

G F m 5 B 

48-7T 3 

4i' 3 (w 

- 1)(™ : 

dr 

dw 

(B- 

-+ Dili/) = 

= |K h | 2 

G F m B 

487T 3 

2r 3 (w 

+ 1 ){w- 

dr 

dw 

l B- 

+ D* 2 lv) = 

= IK b | 2 

G F m B 

487T 3 

2 r 3 (w 

+ l)(wr 


(3.7) 

(3.8) 


2r 3 (w + l)(w 2 — l) 3 / 2 [(w + 1)(1 — r) 2 + 3u>(l + r 2 — 2rw)}\T 3 / 2 (_w)\ 2 (3.11) 


We have also computed dT/dEi where Ei is the energy of the final lepton. We do not think it very useful to 
write down the lengthy resulting formulae, we will show some plots later on. 


4 Bjorken sum rule 

It has been demonstrated in || and J| that the heavy quark models a la Bakamjian-Thomas satisfy exactly the 
Bjorken sum rule i, § thanks to the closure property of the Hilbert space. In this section, we would like to make 
this statement concrete, while providing a flavor of the value and shape of the Tj functions. 

In table [[[ we check the derivative Bjorken sum rule, obtained by a derivation of both sides of full Bjorken 
sum rule at w = 1 : 



We have given the values of Tj(w) for w = 1.0 and w = 1.5 to give some feeling of the decrease of Tj with w. Adding 
only the ground state l = 1 mesons in the r.h.s gives a contribution about 20 % below the l.lr.s, p 2 (35 % in the case 
of the linear + Coulomb potential of VD). We have computed the contributions from the first radial excitation of the 
l = 1 mesons. Now the discrepancy is of a few percent except for the VD, still about 20 % too low. Finally, adding 
all the l = 1 eigentstates in our truncated Hilbert space, we get a 1 % agreement with the directly computed p 2 , thus 
confirming that the models a, la BT do incorporate naturally the Bjorken sum rule. 

The larger value of p 2 in ISGW model is related, via the derivative Bjorken sum rule, to both a larger ti/ 2 (1) 
and a larger t 3 / 2 ( 1). 


The full Bjorken sum rule writes || 


i = ^l£HI 2 + (™-i) 


Tlmax 

E 

. n= 1 


w 2 — 1 




rimax 


(4.2) 


where the dots represent the l > 1 states. In ||, n m ax corresponds to some scale dependent cut-off, which does not 
exist in our case since we did not introduce QCD radiative corrections. We will use for n max the higher state in our 
truncated Hilbert space. Our results are exhibited in tabic |^. 


For w = 1, the sum rule is a triviality from the normalization £(1) = 1 and the vanishing of other contributions. 
The larger w, the more room is left for l > 1 states. However, even for w = 1.5, which is beyond the end of the physical 
domain (w ma x — 1-32 for D**’ s), less than 20 % is due to / > 1 except for the non-relativistic Hamiltonian. Again, as 
for p 2 , we note a general agreement between the three models using the relativistic kinetic energy ®, and a sensibly 
larger (~ 25 %) need of l > 1 states for the ISGW model which uses a non-relativistic kinetic energy. 


5 Numerical predictions for the functions £(w) and Tj(w). 

5.1 The difference 73/2(1) — ti / 2 ( 1 ) 

Table []] clearly shows a dominance of t 3 / 2 (w) over Ti/ 2 {w) (w = 1 — 1.5,) of more than a factor 2, except for the 
non-relativistic ISGW model. This feature is at the origin of a predicted suppression of the semi-leptonic decay into 
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GI 

VD 

CCCN 

ISGW 

7-1/2 (1) 

0.2248(1) 

0.13(1) 

0.059(1) 

0.3424(1) 

7-3/2 (1) 

0.5394(1) 

0.43(2) 

0.515(1) 

0.5918(2) 

71/2(1-5) 

0.1151(1) 

0.10(1) 

0.042(2) 

0.2116(2) 

7 - 3/2 (1-5) 

0.2767(1) 

0.23(1) 

0.270(1) 

0.2700(1) 

ln/ 2 (i)| 2 

0.051 

0.017 

0.004 

0.117 

2|73/ 2 (1)| 2 

0.582 

0.37 

0.531 

0.701 

ground state BSR 

0.882 

0.64 

0.79 

1.068 


0.1965(3) 

0.15(1) 

0.249(2) 

0.2290(2) 


0.2147(1) 

0.25(2) 

0.216(1) 

0.2522(4) 

l T l/2 (1) 1 2 + 2| 7 - 3)2 (1) 1 2 

0.130 

0.15 

0.16 

0.18 

BSR 

1.023 

0.97 

0.97 

1.283 

P 2 

1.023(3) 

0.98(2) 

0.97(2) 

1.283(1) 


Table 1: Results for Tj( 1) and Tj( 1.5) form factors into the D** ’s (lowest 1=1 states). The errors in parenthesis 
is for truncation ( estim ated by varying n max ). The line “ground state BSR” indicates the r.h.s of the derivative 
Bjorken sum rule (4.1) in which the t’s have only been taken for the lowest orbitally excited final states. Adding 
only the contribution of the first radially excited 1=1 final states (indicated by the superindex (1)) almost fills 
the gap with p 2 . The line “BSR” gives the contribution to the r.h.s of the derivative Bjorken sum rule from all 
the states in our diagonalization procedure. The agreement with the l.h.s., p 2 , is almost perfect. 



GI 

VD 

CCCN 

ISGW 

w=1.0 

w=l.l 

w=1.2 

w=1.3 

w=1.4 

w=1.5 

1.000 

0.9864(2) 

0.9546(5) 

0.913(1) 

0.870(2) 

0.825(3) 

1.000 

0.988(2) 

0.962(3) 

0.928(5) 

0.890(6) 

0.85(1) 

1.000 

0.988(2) 

0.959(3) 

0.921(5) 

0.879(2) 

0.83(1) 

1.000 

0.9793(1) 

0.9337(1) 

0.8781(1) 

0.8201(1) 

0.7633(1) 


Table 2: r.h.s of the full Bjorken sum rule (4.2) as a function of w adding all the l = 0 and l = 1 states, 
error is for the truncation. 


The 


D\(l/2)lv and D^lv. It should be noted that this r 3 / 2 (w) i$> Ti/ 2 (w) is only valid for the groundstate l = 1 states. 
Both r’s are of the same order for the first radially excited l = 1 states. 

The difference between t 3 / 2 and t 3 i 2 is not due to the difference between the 1/2 and 3/2 internal wave 
functions at rest coming from the spin-orbit force, which is small and has a rather moderate effect. It is essentially 
due to the relativistic structure of the matrix element in terms of these wave functions. One can gain understanding 
of this effect at least near w = 1 by returning to the analysis made in 0, in a frame where v and v' are small. There, 
the current has been shown (eq.(30) in j^|) to have one spin-independent piece and one which is spin-dependent. 
The latter, corresponding to the light spectator quark Wigner rotations, is responsible for the large difference between 
73 / 2 ( 1 ) and ri/ 2 (l) : 


73/2(1)-7 i /2 (1) 


1 

2V5 


dpp 2 p <j>'(p) 


P 


Po + TO 


Tip) 


(5.1) 


where (f> is <f> 3 / 2 ~ </> 3 / 2 . In a non relativistic expansion, this Wigner rotation effect is small, of relative order 0(v 2 /c 2 ) 
with respect to the main, spin-independent term. That it may be large in spite of this illustrates the fact that the 
system is ultrarelativistic in the GI case, while it is less relativistic in the ISGW case. For p 2 , the Wigner rotation 
effect was found rather small, partly due to a small algebraic coefficient, and partly to the fact that in p 2 it is squared. 





















Figure 1: Plots of £(w). The dots are £(w) calculated for w = 1.0 — 1.5 The curves are functions (2/(1 + w)) 2p f u 
with pj 1t fitted to the above-mentioned points. From upper to lower the curves correspond to VD, CCCN, GI 
and ISGW models. The numbers indicate fjj it for GI, VD, CCCN and ISGW in this order. 


Tau32: sigma32 A 2 - 1.49549, 1.39264, 1.44502, 1,7582 



Taul2: sigmal2 n 2 = 0.831233, 0.57188, 0.733199, 1.07791 



Figure 2: Plots of t 3 / 2 (w) to the left and ^^(w) to the right. The dots are Tj(w) calculated for w = 1.0 — 1.5 
The curves are functions Tj (1)(2/(1 + w)) 2<T i with o 2 fitted to the above-mentioned points. From upper to lower 
the curves correspond to ISGW, GI, CCCN, VD (ISGW, GI, VD, CCCN) for r 3 / 2 (d/i)- The numbers indicate 
a 2 for GI, VD, CCCN and ISGW in this order. 
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5.2 Shape of the functions $,(w) and Tj(w ) 

In order to analyze the shape of the function f(w), we have fitted the predictions of our models to several standard 
one parameter analytic formulae ||17|. 


£nr(w) 

2 

= exp 

-(2 P 2 

U) + 1 


^ISGW-nonrel ) 

= exp [—p 2 (iL 

-1)]. 

£pole(^) 

( 2 \ 2p2fit 

V 1 + w) 

5 

£lin (w) 

= 1 — p 2 {w~ 

!)• 


w — 1 
w + 1 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

{nr(w) has been proposed in (h|, £ po ie(w) in [ftT) and CiSGW-nomei(w) in It should be repeated that fisGW-nonrei(w) 
is totally different from what is referred to as the ISGW model throughout this paper. We only take from the ISGW 
model the mass operator, and compute f(w) consistently in the framework of the BT method. As a result, as we 
shall see, our £(w) deduced from the ISGW model is much better fitted by £ po ie(w) than by £isGW-nonrei(w). We have 
labelled the latter by ISGW-nonrel in order to avoid any confusion. 

We fit p 2 respectively to these four analytic shapes and estimate the quality of the fit by computing: 


5 = 


Y [£K) -£fit{wij? 


1/2 


(5.6) 


^ Wi —1 — 1.5 


where ffu is one of the functions in (5.2)-(5.5) and where uy runs from 1 to 1.5 by steps of 0.1. The results are 

<5nr — (5 to 9) 10 -3 , fecw-nonrei ^ (H to 14) 10" 3 , «5 poIe ~ (2 to 7) 10~ 4 , (S lin ~ (3 to 6) 10- 2 (5.7) 

£ P oie( w ) w i ns by f ar - In fact, in our case, it would deserve the label “dipole” since p 2 ~ 1 and we may claim that the 
BT models with relativistic kinetic energy predict 


- 


(4 m B m D (»)) 2 


{q 2 - (m B + m D (,)) 2 ) 2 


(5.8) 


although, clearly, this singularity does not correspond to an exchanged particle. This dipole behavior is to be compared 
to the pole-type behavior which results, for w —> oo, from a naive Lorentz contraction, see section 6 in |]l9f . Notice that 
the form factors Ai(q 2 ) and fo(q 2 ) have one “pole” less than £(w) (eq (17) in ||l^|), while the other form factors behave 
like £(w). Hence, the result (|5.8| ) favors a pole/dipole behaviour of the form factors, contrary to the constant/pole 
advocated in [[jv) on the basis of the naive Lorentz contraction. In fact other arguments were also used in jf9| and an 
update of the latter discussion would be welcome in view of the theoretical and experimental progress. 

In JTl[ ] we have argued that, when using the relativistic kinetic energy, the ground state wave functions are 
rather close to ground state Coulomb wave function (see eq. (14) anf figure 1 in ]TT||). It happens that in the limite of 
a vanishing light quark mass, the ground state Coulomb wave function gives exactly £(w) = (2/(1 + re)) 2 , i.e. our best 
fit fpoie(w) with p 2 = 1. Clearly the “dipole” behaviour of f(w) may be traced back to the Coulomb-like shape of the 
wave function, which in turn stems from the kinetic energy & It was also noticed in |ll| that this coulombic wave 
function is similar to the NRQCD lattice results. 

The fits with £ po ie(w) are shown in figure 1. The fitted p 2 are presented in table || It appears that they agree 
very well with the exact slope p 2 except for the VD model where a 6% difference is to be noted. We do not know 
whether the latter difference has to do w ith th e str ong si ngu larity of the Coulomb potential in VD model. Needless 
to say, the other fitting functions in eqs. (52), (52) and (|5.5|) le ad to values of p'j it which disagree with p 2 by several 
percent in the case of £nr(w), up to 30% in the case of eq. ( |5.5|) . The latter point has already been discussed in [|TT|1 . 

In a recent paper S. Simula |2C|| has computed the Isgur-Wise function from the light-front constituent model 
of |2l[ based on the GI Hamiltonian. Not only does he get p 2 = 1.03 in the infinite mass limit, but the detailed 
agreement of his solid line in figure 4 with our results for GI model (figure [I]) is very striking. This confirms that the 
models are indeed equivalent as discussed in || . 

Expanding 


f(w) = 1 — p 2 (w — 1) + c(w — l) 2 + ... 


(5.9) 
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GI* 

VD 

CCCN* 

IS GW 

p z fit 

p 2 

n/2(i) 

fJ l/2 

73/2(1) 

2 

fJ 3/2 

1.020 

1.023(3) 

0.2248(1) 

0.83 

0.5394(1) 

1.50 

0.92 

0.98(2) 

0.13(1) 

0.57 

0.43(2) 

1.39 

0.97 

0.97(2) 

0.059(1) 

0.73 

0.515(1) 

1.45 

1.287 

1.283(1) 

0.3424(1) 

1.08 

0.5918(2) 

1.76 


Table 3: Results of the fit of £(u>) by [2/(1 + w)] 2p f it , and of Tj(w) by Tj(1)[2/(1 + w)] 2a J . For comparison 
the exact p 2 is repeated. We have starred our preferred models: GI and CCCN which have a relativistic kinetic 
energy and have a regular potential at r —> 0. 


the different analytic functions in (5 1 2)-(5 1 5) lead to 

p 4 + p 2 - 0.25 


CNR = 

For GI (VD) model it leads to 


p 2 (l + 2p 2 ) 


‘ i ClSGW-nonrel — ^ ’ Cpole — ^ 


p 2 = 1.02, CNR = 0.90, ClSGW-nonrel = 0.52, C po le = 0.77 (GI) 


p 2 = 0.92, CNR = 0.76, ClSGW-nonrel = 0.42, C po le = 0.65 (VD) 


(5.10) 


(5.11) 


In the (p 2 ,c) plane, our best fit (p 2 ,c po i e ) for GI (VD) is slightly above (just on) the ellipse limiting the allowed 
domains shown in figure 1 of |22| . 

In figure 2 we also show the fits of Tj (w) according to 


Tj(w) 




/ 2 \ 2<T **‘ 
\l + w) 


Notice the stable and rather large slopes o 2 ~ 1.5 for the dominant j = 3/2 channel. 


(5.12) 


6 Predicted branching ratios 


6.1 Differential decay widths 


i dr 


our predictions in the GI model for — —— of the semileptonic production of non-orbitally excited 

I aw 


We show in figure | 

D and D* mesons, eqs. (343), (3/7). We compare the latter with DELPHI measurement for B —* D* (fig. 10 in 
multiplied by a factor such that the total number of good events is adjusted to the integral of our predicted partial 
width. The agreement is quite satisfactory. 

Moreover, as promised in a preceding section, we also computed the leptonic spectra dT/ dE\ for each type of 
c-meson produced in the B decays. There is no analytical formula available, as it is the case for the dT /dw widths: 
numerical calculations have to be made. So we are presenting these results in the following plots: the dynamics used 
is, here again, the GI model and we have only considered the non radially-excited D 1 D*, Dg, D\, D\ and mesons 
(though it is possible to compute the same quantities for the radially excited ones). 

In figure [| the leptonic spectrum of the l = 0 D and D* mesons are represented. In figures || and |] the leptonic 
spectra of the l = 1 D** mesons are plotted for the j = 1/2 multiplet and for the j = 3/2 multiplet respectively. We 
can already notice two features that will be discussed in the next section: 

• the j = 1/2 multiplet seems to be outweighed by the j = 3/2 one 

• the contribution of the meson is bigger than the contribution of the D\ meson 
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0.2 



Figure 3: yyg, of the semileptonic production of non-orbitally excited D, eq. (dashed line) and D *, eq. 

( |3. 1\ ) (solid line) mesons, compared with DELPHI measurements for D *, normalized to adjust the total number 
oi experimental events to the integral of our theoretical curve. The experimental numbers correspond to bins in 
w of size 0.1. 




Leptonic spectrum B —> D I v Leptonic spectrum B —> D’ I v 


Figure 4: Leptonic spectra - of the semileptonic production of l = 0 D and D* mesons 

T dEi 
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Figure 5: Leptonic spectra L/tt of the semileptonic production of the first j = 1/2 multiplet of D* 


mesons. 




Figure 6: Leptonic spectra A -/£- of the semileptonic production of the j = 3/2 multiplet of D** mesons. 







6.2 Comparison with experiment 


The summary is presented in table ||. The agreement is quite satisfactory for D^*\ 

Concerning the orbitally excited states, some work has to be done to compare the different experiments between 
themselves and with our models. Experimental numbers are typically provided as branching ratios Br{b —> Dj {—> 
D* + n~)lv) (Dj being any of the two observed orbitally excited D mesons: £>2(2460) or Z?i(2420)). To reduce this to 
the branching ratio B —> Djln, we use the nominal production ratio of b jet: B + : B° : B s : 5-baryon = 0.4, 0.4, 0.12, 
0.08 @. 

The branching ratios into D*tt are very different for and D\. Indeed, the decay D\ —> Dn is forbidden 
by angular momentum and parity conservation. We have assumed that these two channels D^*’n saturate the Dj 
decay. This seems a plausible assumption, although it has been argued Q that three-body, p and/or p decays might 
be important. This means that we assume Br(D\ —► D*n) = 1, and from isospin Clebsch-Gordan coefficients, the 
decay into charged pions Br(D± —> D* 7r ± ) = 2/3. Concerning the -DjJ decay, we have two pieces of information: 
one from experiment, |j(J indicates that .Br/DJ; (2460) —> D + 7r _ )/Br (£>2(2460) —> D* + n~) = 2.2 ± 0.7 ± 0.6, the 
other one form theory, the value 2.3 predicted by HQET |33|]. Since both agree we take this ratio to be 2.3 implying 
Br(D 2 —> Dn) = 0.7 and Br(D 2 —> D*n) = 0.3, with an additional C.G. coefficient of 2/3 for final charged pions. 

We have treated the above mentioned estimates for the production ratios b —> B and for the branching ratios 
Dj —> D*n as if they were exactly known. We do not feel able to discuss in a reliable way the error attached to them. 
As a conclusion the experimental errors mentioned in table || do not incorporate this uncertainty. They underestimate 
the error. Of course, every time the experimental groups have already performed one of the above mentioned steps, 
we use their estimate. Then the quoted error is more realistic. For example CLEO directly provides the B~ branching 
ratios. 


Looking at the table |s|, it appears that the experimental numbers are still rather scattered. The difference in 
OPAL results between charged and neutral D\ is quite surprising. Our models agree with the smallest one, D 2 °. the 
other experiments give only upper bounds for D 2 . None is in clear disagreement with our models. It is standard to 
look for the orbitally excited D mesons via their D*n decay. The lower branching ratio of D\ —> D*n reduces the 
number of observed D 2 compared to £>i, allowing only upper bounds, which are not the sign of a small Br(B —> D 2 lv) 
but only of a small Br(D 2 —> D*n). This has to be kept in mind. 

The results for £>i are rather scattered: ALEPH and CLEO find small values, DELPHI and OPAL larger ones. 
However the discrepancy is less than two sigmas. We agree with the small values, ALEPH and CLEO. 

A special comment is due concerning the ratio Br(B —* D 2 )/Br(B —> D{). In the heavy mass limit, the 
two decays depend on one Heavy-quark universal form factor: t 3 / 2 (' w). Of course, since the kinematic factors differ 
between B —> D 2 and B —> Di (compare eqs. (3.11) and ( }3.1C| )), the ratio of the partial widths slightly depends on 
the detailed shape of the function t^/ 2 {w) which is model dependent. However, this dependence is not very important 
as can be seen from the stability of the ratio in table [|. To get some feeling of this dependence we have considered two 
extreme shapes for r 3 / 2 (ru): a constant and a very steep decrease ( 0 - 3/2 = 3). The ratio Br(B —> D%)IBr(B —> D\) 
ranges from 1.43 to 1.7. We may thus consider a ratio of ~ 1.55±0.15 as a prediction of HQS. Only DELPHI provides 
this ratio, which agrees with HQS but within large errors. 


Channel 

Gl* 

VD 

CCCN* 

IS GW 

exp. 

B -> Din 

B -> D*lv 

B -► D\lv 

B -> D\(%)lv 
B -► DlU)lu 
B^D* 0 (l)lv 

B—DHf)/!/ 

P 2 

2.36 

6.86 

0.70 

0.45 

0.07 

0.06 

1.55 

1.023(3) 

2.51 

7.19 

0.47 

0.29 

0.02 

0.02 

1.62 

0.98(1) 

2.45 

7.02 

0.65 

0.42 

0.004 

0.004 

1.55 

0.97(1) 

1.94 

6.07 

0.77 

0.49 

0.13 

0.11 

1.57 

1.283(1) 

1.6 ± 0.8 
5.3 ± 0.8 


Table 4: Semileptonic branching ratios, in percent, taking a B life-time of 1.62 10~ 12 s. The last line is the 
prediction for p 2 . The column on experiment refers to the rather established and stabilized results reported in 
the PDG table. For more recent experimental results, see the following discussion and table 5. We have starred 
our preferred models as in table 3. 
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Channel 

Theory 

PDG 123D 

ALEPH 125(1 126| 

DELPHI lp7|| ]2EJ 

OPAL 1 291 

CLEO 1311 

B -» Dlv 

B -> D*lv 

B -> Dfh/ 

B -* D* + h/ 
B-*Di(%)lv 
B -► Dl\l)lv 
B^D* 0 (l)h/ 

b^TdrI )TT 

1.9 - 2.5 

6.0 - 7.2 
0.45 - 0.80 

0.45 - 0.80 
0.30 - 0.50 
0.0 - 0.07 
0.0 - 0.06 

1.55 - 1.62 

1.6 ±0.7 
5.3 ±0.8 

2.35 ±0.20 ±0.44 

5.53 ±0.26 ±0.52 

< 0.94 

< 0.54 

0.74 ±0.16 

5.47 ±0.16 ±0.67 
< 6.25 

1.5 ±0.55 

1.55 ± 1.1 

0.88 ±0.35 ±0.17 

2.25 ±0.65 ±0.51 
2.0 ±0.6 ±0.5 

< .85 

.49 ± .14 ± .06 


Table 5: The first column gives the range of our theoretical predictions. We have assumed a Br(b —> B + ) = 
Br(b —> B° ) = 0.4 in order to compare the experimental numbers referring to Br(b —> Dj..) with our prediction 
concerning Br(B —> Djlu). Whenever experiment gives a branching ratio —> DjlX we have assumed X to 
he always a neutrino. We have furthermore assumed Br{D\ —> D*n) = 3.33 Br(D 2 —> D*tt) from the ratio 
BrfD^ —*• D-k)IB r(T >2 —> .D*7r) ~ 2.3 and neglecting other D\ decay channels. Unable to make a reliable 
estimate of errors on all these estimates about Br(b —> B) as well as about the Dj decay branching ratios, we 
have chosen to take them with zero error. It results that the experimental errors on the numbers given in the 
table are indeed larger than indicated. The upper bounds are at 95 % CL, except for the caveat in the previous 
sentence. 


However, if we take strictly this HQS prediction Br(B —> D^lv) / BrfB —> Dili/) > 1.4, it appears that the 
central value for Br(B —> Dili/) in ALEPH and OPAL leads to a prediction for Br(B —> D^lv) which overcomes the 
upper bound of the same experiment. In view of the general uncertainty on these numbers, we can only take this as 
a faint indication that some of our hypotheses might need to be reconsidered. 

Let us first consider the effect of neglecting all decay channels of the D **'s except into D^n. Let us assume 
that some other channels contribute with a partial width to both D% and D\ decays, assumed to be equal to the partial 
width into D*n times r. As a result we should divide in table [| the experimental numbers concerning Di by 1 + r 
while those concerning D | will be divided only by 1 + 0.3 r. it results that this correction would bring the numbers in 
better agreement with the HQS prediction Br{B —* D^lv) / Br(B —> Dili/) > 1.4. At this stage we conclude that it 
would be prematurate to claim that experiment contradicts this HQS ratio. 

Let us still assume that there is some indication of a needed correction to the strict HQS relation Br{B —> 
D^lv)IBr(B —> Dili/) > 1.4. This leads to consider 0(l/m c ) corrections. A mixing of the Di with the j = 1/2 meson 
is quite plausible, but it would worsen our prediction since, due to the very small decay amplitudes into the j = 1/2 
states, it would lessen our prediction for Br(B —■> Dili/). We are thus lead to consider direct 0(1/m c ) corrections in 
the decay amplitudes, which have indeed been found to be dramatically large in Q. The authors of the latter paper 
obtain Br(B —> D%lv)/Br(B —> Dili/) ~ 0.8 for constant Tj{w)’ s. 

Unhappily, no direct check of the predicted suppression of the semileptonic decay into j = 1/2 or bitally excited 
D's is presently available. 


6.3 Comparison with other models 


Quark models. Among the calculations of b —> c exclusive semi-leptonic decays which have been performed in the 
framework of quark models, one must distinguish between: 1) Calculations in frameworks equivalent to the present 
B-T formalism |36), Q, Q. 2) Calculations with direct identification of some scaling functions to NR overlaps : 
j35|. [|8|. 3) Calculations in a family of models with approximate Lorentz boost |3!|, which are close in spirit to 
our old quark model Jl9| . 4) Calculations which do not deduce the q 2 dependence of form factors from wave functions, 
but from some simple or multiple pole Ansatz ; pi) (these models have been discussed in Q), C3] ; or still some 


other Ansatz 116 


1) Calculations in frameworks equivalent to the present B-T formalism. In spite of the fact that the authors 
do not give general formulae, we have found that the results of ref. |3(J concide with the ones of the B-T method for 
the particular case they consider (tuq = oo limit, choice of a particular frame, expansion around w = 1, harmonic 
oscillator wave functions). Therefore we state that they have obtained by intuitive arguments essentially the B-T 
formulation for mQ = oo. The quite sizeable differences observed with respect to our present quantitative predictions 
for p 2 and ti/ 2 ( 1 ), 13 / 2 ( 1 ) come from our use of more realistic wave functions: the full numerical solution of relativistic 
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spectroscopic models (using the KE of eq (2.1)), instead of the non relativistic ISGW spectroscopic model with 
harmonic oscillator (h.o.) approximation, yield lower results for p 2 . 

As to the work of ^0|, |l7|, it is a full calculation of the B —> D^iv form factors with the GI spectroscopic 
model, in a null-plane approach, which, as we have checked (ref. [Q ), is the application of the B-T method in the 
particular P —> oo frame. As we have suggested [Q, because of the demonstrated covariance of the niQ = oo limit, 
the results should then coincide in this limit with our own, provided the same spectroscopic model is used. Indeed 
our £(u>) agrees within the computational uncertainty with the prediction by these authors when using the same GI 
model. 

2) Calculations with direct identification of some scaling functions to NR overlaps. In J35|, J38|, a simple 
intuitive recipe has been used, consisting in the direct identification of the scaling functions of Isgur and Wise with 
usual non-relativistic form factors. It results in quite different expressions from what we advocate througout this 
paper: deducing the scaling functions from their relation to the current matrix-elements, the latter being calculated 
in a relativistic approach, as is proposed by the B-T approach and the one below (point 3)). 

3) Calculations in a family of models with approximate Lorentz boosts. In J39| , |4(]] , the approach is essentially 
the one of our old model developped in the 70’s (in the h.o. case, one ends with exactly the same £ that we had found; for 
a recent reference, see |l^|). It consists in calculating the current matrix elements between wave functions for states in 
motion obtained through a Lorentz transformation of rest-frame wave functions, including the Lorentz transformation 
of coordinates and the standard transformation of spinors in Dirac representation. The matrix elements are calculated 
in the equal-velocity frame where great simplifications occur. To get the equal time wave function from the one at 
rest, one has to assume that the dependence on the relative time is weak in this latter frame (this corresponds to the 
assumption of factorisation of the wave function in light and heavy degrees of freedom in (39|, [Q. One then ends 
with what we have called the “naive Lorentz contraction” prescription for spatial wave functions. This prescription 
can be seen to be too inaccurate to calculate form factors. In the B-T formalism, there appears a definite and more 
consistent Lorentz transformation of momenta. The “naive Lorentz contraction” prescription consists in peforming 
in the result of the BT boosts a linear approximation in the dependence on the light quark spatial momenta p 2 in 
the hadron rest frame . The resulting difference can be seen in the behavior of the function f,(w). With the simpler 
Lorentz contraction prescription, t;(w) = 2I(w)/{\ + w) where the overlap factor I(w) —> constant when w —> 
Hence f(w) ~ (2/ (w + 1)) for large w. In the B-T formalism, as we have argued after eq. (0U(«O - (2/ (w + l)) 2p , 
with a power of dynamical origin. 

Another point is deserving mention as regards models |]39| , j|4(i|| . Within a static approximation where the light 
quark wave functions at rest would have no small Dirac components (’’non relativistic” case of JfJ), and assuming, as 
is found to a good accuracy, that the spin-orbit effect is small, one gets Ti/ 2 { 1) — r 3 / 2 (1). On the other hand, if one 
uses wave functions with small components, like for instance solutions of a Dirac equation (”DESC” of jif|), one has 
r i/ 2 (l) 7 ^ t 3 / 2 (1), and indeed, one finds, similarly to our B-T result, T( 1/2) << T(3/2). In fact the B-T formalism 
corresponds to having small components given by the free-quark Dirac spinor structure. 

Apart from particular features, one must emphasize the theoretical advantage of the B-T approach that 
comes from the demonstration of exact general properties such as unitarity of the transformation from rest 
frame to states in motion, or, in the toq = oo limit, Isgur-Wise scaling, full covariance and duality. In the 
other approaches, either they do not hold, or they are enforced by formal procedures, for example formal covariance 
can be obtained by identifying invariant form factors with form factors in one particular frame like the equal-velocity 
frame, or they are only approximate. 

Finally, one must still mention the following important point: the above discussion has been maintained within 
the toq = oo limit. For our B-T approach, as well as for the others, it is however quite possible to treat the finite 
mass case, and it is indeed one interesting feature of quark models that they yield definite 1 /toq effects. We have 
only refrained to do so because our approach loses several of its nice properties, in particular covariance. Detailed 
predictions for finite mg are given in the ISGW2 model j42|. This model falls in the category 4), in that the form 
factors do not really derive from a full quark model calculation. One interesting finding of the ISGW2 calculations 
||42f is that they predict T(2 + ) < T(l^ 2 ), in contradiction with HQS. This may be due to the finite mass treatment 
and would be in agreement with the conclusions of |Q. It may also be mentioned that for the ISGW1 model ( ]lGtl ), 
r 3 / 2 (l) is sensibly smaller than ours (0.3 against 0.55). 

QCD fundamental methods. By fundamental methods, wee mean lattice QCD and QCD sum rules. As to 
the comparison of their results with ours as regards the elastic Isgur-Wise function the discussion has been made 
in [|ll], and the convergence is encouraging, even as concerns the sensitive p 2 parameter. For orbitally excited states, 
only QCD sum rules give predictions ( |43| ), and they lie well below the ones of our quark models for ununderstood 

5 In the harmonic oscillator case £{w) = £nr(w) (see eq. (|~>.i])) 
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reasons (with toq = oo, Br(B —> D\[3/2) + lv) = 0.1% versus 0.3 — 0.5% in table ||). 


7 Conclusion 

We have applied the Bakamjian-Thomas (BT) method in the infinite mass limit to meson wave functions fitted to the 
experimental meson spectrum by several groups: GI, VD and CCCN, 0-0, using a relativistic kinetic term and 
ISGW 0] a non relativistic one. Whichever set of wave function is used, the BT method insures that the models are 
covariant and satisfy Heavy Quark Symmetry i-0- 

We have computed the invariant form factors £(w) and Tj(w). The Bjorken sum rule, which has been proven 
to be valid in BT models ||, |(J, has been checked practically in this paper. The derivative Bjoren sum rule for w = 1 
is almost saturated by the ground state 1 = 1 and the first l = 1 radial excitation (the latter contributing for ~ 15%), 
except for the VD model, in which ~ 20% is contributed by the n > 1,1 = 1 excitations. The full Bjorken sum rule 
has also been checked for 1 < w < 1.5. The missing part, corresponding to the l > 1 states, increases with w up to 
15% for w = 1.5 (25% for the ISGW model). 

The slope p * 2 is 1. ± 0.05 when relativistic kinetic energy is used, ~ 1.3 when the non-relativistic kinetic energy 
is used_0|. The £,(w) function is very well fitted by the function (2/(w + l)) 2p which means a dipole like behaviour, 
eq. ( |5.8| ). The latter result may be traced back to the Coulomb-like shape of the ground state wave function when 
using a relativistic kinetic energy 0. The functions Tj are dominated by t 3 / 2 , with r 3 / 2 (l) — 0.4 — 0.6. The latter 
decrease like (2 /{w + l)) 3 4 5 . 

We have computed, for the decays B —■> Dlv 1 D*lv and the four B —> D**lu , the differential decay widths as 
a function of w, figure |], and also as a function of the charged lepton energy. Figures [i||(] show the latter. The decay 
widths to j = 3/2 are about one order of magnitude smaller than the ones into 1 = 0 states, and the ones into j = 1/2 
are still one order of magnitude below. 

We have also computed the semileptonic branching ratios. The ones into D and D* are in good agreement 
with experiment. For the D** experimental data are still rather scattered. Our models predict tiny branching ratios 
into j = 1/2 states, but this cannot be checked yet. For sure these decays have not been seen yet, but their absence 
may be explained by their broad widths, without invoking a suppressed branching ratio. Our predictions for D 2 are 
below the experimental upper bounds, although close to them. Our predictions for the Di(j = 3/2) do not disagree 
with experiment taking into account the wide spreading and large uncertainty of experimental results. 

Our models being all within the strict infinite mass limit, we predict the relation Br(B —> D 2 )/Br(B —> 
D\(j = 3/2) = 1.58 ± 0.05 which is mainly a consequence of HQS. Indeed, going beyond the models studied in this 
paper, if we consider arbitrary but “reasonable” shapes of t 3 / 2 (w ) we find that HQS predicts Br(B —> D^)/Br(B —> 
Di(j = 3/2) = 1.55±0.15. When comparing the experimental figures for these decays, there might be some indication 
that this ratio is indeed lower. This could be an indication of a large l/m c correction which enhances the D\. 
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